KUNS 2288 



Catastrophic Instability of Small Lovelock Black Holes 

Tomohiro Takahashi and Jiro Soda 
Department of Physics, Kyoto University, Kyoto, 606-8502, Japan 

(Dated: August 19, 2010) 

Abstract 

We study the stability of static black holes in Lovelock theory which is a natural higher dimen- 
sional generalization of Einstein theory. We show that Lovelock black holes are stable under vector 
perturbations in all dimensions. However, we prove that small Lovelock black holes are unstable 
under tensor perturbations in even-dimensions and under scalar perturbations in odd-dimensions. 
Therefore, we can conclude that small Lovelock black holes are unstable in any dimensions. The 
instability is stronger on small scales and hence catastrophic in the sense that there is no smooth 
descendant. 

PACS numbers: 98.80. Cq, 98.80.Hw 
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I. INTRODUCTION 



The possibility of higher dimensional black hole creation at the LHC would be the most 
fascinating prediction of the braneworld with large extra-dimensions Nowadays, it is 
fashionable to study higher dimensional black holes. The study of higher dimensional black 
holes is not a straightforward extension of the 4-dimensional one. In fact, in higher dimen- 
sions, there are many black objects with non-trivial topologies such as black rings. This 
seems to require the stability analysis of higher dimensional black objects. Moreover, in 
higher dimensions, Einstein theory is not the most general theory of gravity which contains 
terms only up to the second order derivatives in the equations of motion. The most general 
theory in this sense is Lovelock theory Q . Thus, it is worth extending the stability analysis 
to this general Lovelock theory. 

The stability of higher dimensional black holes has been intensively studied since the 
seminal papers by Kodama and Ishibashi 3j. It is important to study various black holes 
in Einstein theory because black holes produced at the LHC are expected to be charged or 
rotating. A numerical study of charged black holes has been done Q. Ib investigate the 
stability of rotating black holes, a group theoretical method is developed [5[. The method 



is used to study the stability of squashed black holes js-8| and 5-dimensional rotating black 
holes with equal angular momenta {9]. The stability of a special class of rotating black 



holes in more than 5-dimensions is also studied 
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10Hl2|. Recent development of numerical 



stability analysis is remarkable 13l-ll7j]. It has been shown that rapidly rotating black holes 
are unstable. 

As we mentioned already, in higher dimensions, we need to extend Einstein theory to 
Lovelock theory. In addition to the theoretical requirement, we have a physical motivation 
to consider higher derivative theory of gravity. In fact, at the energy scale of black hole 
production, Einstein theory is not reliable any more. It is believed that string theory which 
can be consistently formulated only in 10-dimensions is the most promising candidate of 
the unified theory. We should recall that string theory predicts Einstein theory only in the 
low energy limit 18]. In string theory, there are higher curvature corrections in addition to 



the Einstein-Hilbert term 



. Thus, it is natural to extend gravitational theory into those 



with higher power of curvature in higher dimensions. It is Lovelock theory that belongs 



to such class of theories 



In Lovelock theory, it is known that there exist static 
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spherical symmetric black hole solutions 20[ (and topological black hole solutions are also 



found in 21]). Hence, it is natural to suppose black holes produced at the LHC are of this 
type 22J. Thus, it is important to study the stability of these Lovelock black holes. 

In the case of the second order Lovelock theory, the so-called Einstein-Gauss-Bonnet 
theory, the stability analysis under tensor perturbations has been performed 23j (see also 



an earlier work 
bations 



241] ). The analysis has been also extended to the scalar and vector pertur- 



251 ] . It is shown that there exists the scalar mode instability in 5-dimensions, the 



tensor mode instability in 6-dimensions, and no instability in other dimensions. Although 
Einstein-Gauss-Bonnet theory is the most general theory in 5 and 6-dimensions, it is not so 
in more than 6-dimensions. For example, when we consider 10-dimensional black holes, we 
need to incorporate the fourth order Lovelock term. Indeed, when we consider black holes 
at the LHC, it is important to consider these higher order Lovelock terms 26|. Hence, the 



purpose of this paper is to study the stability of black holes in any order Lovelock theory, 
namely, in any dimensions. We have already shown that Lovelock black holes are unstable in 
even- dimensions under tensor perturbations [27]]. In this paper, we extend previous results 
to vector and scalar perturbations using the master equations we have obtained recently 28]. 

The organization of this paper is as follows. In section [Til we review Lovelock theory 
and explain a graphical method for constructing Lovelock black hole solutions. In section 
IIH] we consider tensor perturbations and show the instability of small Lovelock black holes 
in even- dimensions. This is a review of our previous paper 23] ■ In section ITVl we show 
that black holes are stable under vector perturbations. In section [V] we examine scalar 
perturbations and show that there exists the instability in odd-dimensions if black holes are 
sufficiently small. In section I VI] we present a detailed analysis for Einstein-Gauss-Bonnet 
theory to illustrate our statements. The final section IVHI is devoted to the conclusion. 



II. LOVELOCK BLACK HOLES 



In this section, we review Lovelock theory and introduce a graphical method to reveal 
the nature of asymptotically flat black hole solutions. 

The most general divergence free symmetric tensor constructed out of a metric and its first 
and second derivatives has been obtained by Lovelock {2]. The corresponding Lagrangian 
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can be constructed from m-th order Lovelock terms 



_ fi^l (Tl • • • A m CT m J£ 



(1) 



where R\ a pK is the Riemann tensor in D-dimensions and 5^ S -Z 1 '"^ m ^ m is the generalized totally 
antisymmetric Kronecker delta. By construction, the Lovelock terms vanish for 2m > D. It 
is also known that the Lovelock term with 2m = D is a topological term. Thus, Lovelock 
Lagrangian in D-dimensions is defined by 



(2) 



m=0 



where we defined the maximum order k = [(D — l)/2] and c m are arbitrary constants. Here, 
[z] represents the maximum integer satisfying [z] < z. Hereafter, we set Co = — 2A, c\ — 1 
and c m = a m /m (m > 2) for convenience. Taking variation of the Lagrangian with respect 
to the metric, we can derive Lovelock equation 



1 a 
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(3) 



As is shown in 
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m=l 



there exist static exact solutions of Lovelock equation. Let us 



consider the following metric 



dv 2 

ds 2 = -f(r)dt 2 + — — + r 2 ^ ij dx i dx j 
f( r ) 



(4) 



where is the metric of n = D — 2-dimensional constant curvature space with a curvature 
k—1,0 or -1. Using this metric ansatz, we can calculate Riemann tensor components as 



R 



ti- 



tr 



f 

~~ 5 Rti Rri 



f . 

—Si , Rij 



kl 



(5i h 5j l - S/S/) 



Substituting (jSJ) into and defining a new variable ^(r) by 

/(r) = k — r 2 ip(r) , 



we obtain an algebraic equation 



W\ 



£ 

m=2 



2m-2 



rn 



[(n-p)U r 



p =i 



2A 



nm + 1 



r.n+1 



(5) 



(6) 



(7) 



In (J7J), we used n = D — 2 and /i is a constant of integration which is related to the ADM 
mass as 



M 



2/i7r(" +1 )/ 2 
r((n+l)/2) ' 
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(8) 
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FIG. 1: The intersection between the solid curve and the thin horizontal line determines the 
solution ip = ip(r) for the case n = 4. Apparently, the infinity r = oo corresponds to if) = 0. The 
intersection between solid and dashed curve gives a horizon ru- 



where we used a unit 167rG = 1. 

From (JTJ), it is easy to see that the solution f(r) has many branches. In this paper, we 
want to concentrate on asymptotically flat spherically symmetric, i.e. A = and k = 1, 
solutions with a positive ADM mass \i > because such black holes could be created at the 
LHC. We also assume that Lovelock coefficients satisfy 



a m > 



(9) 



for simplicity. For example, in the case of n = 3, the theory is reduced to Einstein-Gauss- 
Bonnet theory. In this case, Eq.Q reads 



a 2 i> 2 + i> = -^ 



(10) 



which can be explicitly solved as 



-1 ± 



2a 5 



The upper branch leads to the asymptotically flat solution 

.2 



f(r) 



1 + 



r 
2^ 



1 



4a 2y u 



;i2) 
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In the case of n = 4, Eq.flTJ) is reduced to 



3a 2 ip 2 + ip = — ■ 



(13) 



Although it is easy to solve this equation analytically, in FigfTJ a graphical method is also 
explained for this case. However, in the case of n — 5, Eq.flTJ) becomes the third order 
algebraic equation 



We have a formula for solutions of Eq. ()14p . however, the roots are complicated in general. 
Hence, we use a graphical method illustrated in Figj2j Because of the conditions (jHJ), the 
function is monotonic for positive ip. From flTJ), we see the root behaves ip ~ fi/r n+1 or 
f(r) ~ 1 — \ijr n ~ x as r — > oo. Thus, the asymptotically flat solutions belong to the branch 
where ip is always positive. 



FIG. 2: We illustrate a graphical method for n = 5 case. In this case, the third order Lovelock 
theory is most general. Therefore, W[ip] in ([T]) reads a cubic polynomial. In this figure, three roots 
are depicted. Among these roots, only ip > one corresponds to an asymptotically flat solution. 

Now, let us look for the horizon of black holes. The horizon radius of the asymptotically 
flat solution is characterized by f{rjj) = 0. From fl6j), we have a relation i/jh = ip{rn) — l/ r lr 
Using this relation and flTJ), we obtain an algebraic equation 



8a 3 ^ 3 + 6a 2 ^ 2 + rj) = . 



(14) 




asymptotically AdS branch 



asymptotically flat branch 



(n+l)/2 



(15) 



H 
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This determines ipH and hence th- For example, in the case of n = 3, we can easily solve 
this polynomial equation ffT5]) as 

1 



(16) 



/i - 0,2 

Note that if /i < a2, there appears a naked singularity. So, we consider the range fi > 0,2- 
In the case of n=4, since it is a bit complicated to solve Eq.( fl5l) analytically, we present 
a graphically method in FigJTJ Similarly, in FigEl we present a graphical method to solve 
Eq.f ll5p for n = 5. From FigJHJ it is obvious that the range 00 > r > r# corresponds to 
< ip < ipH when f(r) describes an asymptotically flat solution. It is also apparent that 
ipH becomes larger as /1 becomes smaller. 

Remarkably, the nature of black holes depends on the dimensions. In even- dimensions 
n = 2k, dividing (fl5|) by ipn, we have 



£ 

m=2 



2m-2 



p=l 



m—1 
H 



+ 1 = 0. 



(17) 



Near ipn = 0, the left hand side is positive; however, when tpn is sufficiently large, it is 
negative. Then, there exists a positive root somewhere between. This root moves from to 
00 as /i moves from 00 to 0. Thus, there is no restriction for fi in this case. On the other 
hand, in odd-dimensions n — 2k — 1, after dividing by ipn, Eq. (fT5l) becomes 




2fc-2 



k-1 



H(n-p) Ur + E 



P =i 



m=2 




2m-2 



[(n-p) US 



m— 1 



p=l 



1 = 0. 



is) 



In order to have a positive root, we need /i > ^ jll^ta 2 ^ — P) j ■ Hence, we have the lower 
bound for the mass in odd-dimensions. In fact, there exists a positive root in this case, 
because the left hand side of Eq. (fT8|) is positive near ipjj = and negative for sufficiently 
large ipn- Furthermore, It is not difficult to see that this root approaches as fi — > 00 
and approaches 00 as /i -> y jE[p!=L 2 ( n — p)j- Therefore, the root ij)H moves in the range 
< ipH < 00 as p, moves in the range ^ |ripii 2 ( n — P)| < M < 00. 

Finally, we examine the singularity in the solutions. Using the metric ansatz (jl]), the 
Kretschmann scalar R flu \ p R' J ' l/Xp is 



R 



f 



'2 



r ^p = f + + 2n{n - 1) 
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FIG. 3: For n = 5 case, this figure explains a method for calculating ijjjj or ru graphically. The 
positive root of Eq. (115j) ipH can be obtained from the intersection of the solid and dashed curves. 
Since the intersection between the horizontal line and the solid curve gives a solution if) = ip(r), the 
horizon radius ru is determined from the intersection of the horizontal line, the solid and dashed 
curves. 



Thus, this solution has curvature singularities at r = or at the point where derivatives 
of f(r) diverges. For example, / diverges at the point where ip diverges because of /' = 
— 2rvp + r 2 ip' . Taking a derivative of (J7J) with respect to r, we have a relation 

, = (n + l)/i 

W ' r n+2 dipW[ilj} ' 

Hence, a curvature singularity appears at the point d^W = 0. However, for the asymp- 
totically flat branch, W[ip] is monotonically increasing function of ip. Hence, we conclude 
9^VF[^] > 0. Similarly, for asymptotically flat spherical symmetric solutions, we can prove 
that the conditions |1 — f\ < oo, \f'\ < oo and \f"\ < oo are satisfied except for r = 0. 
Therefore, there is a curvature singularity only at r = in the cases we are considering. 
And asymptotically flat solutions have the horizon if the parameter \i is sufficiently large. 
So, these solutions do not have a naked singularity and describe black holes with a mass M 
defined by Eq.©. 
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III. STABILITY ANALYSIS FOR TENSOR PERTURBATIONS 



In this section, we examine the stability under tensor perturbations, which we have 



already studied in 



27] 



We start from the master equation for tensor perturbations 
p2 „ / 2 T" , 2p , A , , (2K + lt )fT" 



~f 2 X ~ UV + — + ff * + > > X = ^, (19) 
y 1 r J (n — 2)r 1 

where a; is a frequency and we have defined a key function 

T(r) = r n -%W[i/)\ . (20) 



In ( 1191) , x is t ne master variable and 74 is eigenvalue of tensor harmonics which is given by 
7t = + — 1) — 2, = 2, 3, 4 • • • ) for k = 1 and positive real numbers for n = —1,0. 

Here, we should recall our assumptions. We assumed the conditions (jUJ) are satisfied. And 
we also assumed spherical symmetry and positivity of the mass, i.e., k = 1 and fi > 0. Then, 
there exists an asymptotically flat spherical symmetric branch which we have considered in 
section HIl Note that T(r) which is defined above is always positive in this branch. 



A. Criterion for Stability 

We will present the condition for the stability of the solutions we are considering. 
As we will soon see, the master equation (1T91) can be transformed into the Schrodinger 
form. To do this, we have to impose the condition 

T'(r) > , (for r > r H ) ■ (21) 

In fact, this is necessary for the linear analysis to be applicable. In the case that there exists 
r such that T' (tq) = and r > r#, we encounter a singularity. Using approximations 
T (r) ~ T"(r )(r — r ) = T" (r )y, f(r) = f(r ) and r = r , (|T9|) approximately becomes 

This shows that near r = vq, \ behaves as \ ~ c\ + C2 logy , where c\ and C2 are constants of 
integration. Hence, the solution is singular at y — for generic perturl 



Dations. The similar 



situation occurs even in cosmology with higher derivative terms 30 



3l(]. In those cases, 
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this kind of singularity alludes to ghosts. Indeed, if there is a region T (r) < outside the 
horizon, the kinetic term of perturbations has a wrong sign. Hereafter, we call this the ghost 
instability. 



When the condition ( 12TI) is fulfilled, introducing a new variable \l/(r) = x{ r ) r \/T' (r) and 
switching to the coordinate r*, defined by dr*/dr = 1/f, we can rewrite Eq. ffl9|) as 

d 2 <5> 



dr* 2 

where 



V t (r(r*))V = ou 2 ty , (23) 



is an effective potential. 

For discussing the stability, the " S-deformation" approach is useful |3|, |23|. Let us define 
the operator 

Hs -^ + v ' < 25 > 

acting on smooth functions defined on / = (r* H , oo). Then, fl23|) is the eigenequation and u 2 
is eigenvalue of Ti. We also define the inner products as 



(Vi.^a) = J <p\<p%&r* • (26) 
In this case, for any ip, we can find a smooth function S such that 

(tp,H<p) = J(\D<p\ 2 + V\<p\ 2 )dr\ (27) 

where we have defined 

D = ^ + S, V = V t + ff-S 2 . (28) 
dr* dr 



Following [23[, we choose S to be 



S = -f^-\n(rVT) . (29) 



Then, we obtain the formula 



(V, = / + (2k + 7 t) jf (j!^)/^ ^ • ( 30 ) 
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Here, the point is that the second term in ( 130]) includes a factor 2k + 74 > 0, but T does 
not include 74. Hence, by taking a sufficiently large 2k + 74, we can always make the second 
term dominant. 

Now, let us show that the sign of dlaT /dr determines the stability. If dlnT /dr > 
on /, the solution (j3J) is stable. This can be understood as follows. Note that 2k + j t > 0, 
then we have V > for this case. That means {tp,Htp) > for arbitrary tp if dhxT' /dr > 
on /. We choose, for example, tp as the lowest eigenstate, then we can conclude that the 
lowest eigenvalue u 2 is positive. Thus, we proved the stability. The other way around, if 
d\nT jdr < at some point in J, the solution is unstable. To prove this, the inequality 

(31) 

is useful. This inequality is correct for arbitrary tp. If d In T' jdr < at some point in /, we 
can find tp such that 

M 2 dlnT' , 

, : : dr < . 32 

, (n-2)r dr v 1 

In this case, ( 130]) is negative for sufficiently large Ik + 7 t . Then, the inequality (13TT) implies 
< and the solution has unstable modes. Thus, we can conclude that the solution is 
stable if and only if dlnT' /dr > on I. 

From the above logic, if dlnT' /dr has a negative region, negative u 2 states exist. There- 
fore, this instability is dynamical. Then, we call this as dynamical instability in order to 
distinguish this from the ghost instability which is caused by negativity of T (r). 

We want to summarize this subsection. If T has negative region outside the horizon 
r > rjj, Lovelock black holes have the ghost instability. Even if T is always positive, 
Lovelock black holes have the dynamical instability if T" has a negative region outside the 
horizon. Therefore, Lovelock black holes are stable under tensor perturbations if and only 
if T' and T" are always positive outside the horizon. 

B. Instability of Small Lovelock Black Holes in Even-dimensions 

In this subsection, we check the sign of T and T" for asymptotic flat solutions f(r). In 
order to see the sign of these functions, it is useful to express them as functions of ip instead 
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of r. Using Eq.flTJ) and its derivative, we obtain 



(fyWfy))^ = -(n + 1)- 



■n+2 



-(n+1 



W[ip] 



(33) 



The above formula can be used to eliminate tp' in T' . The result reads 



T(r) 



„n-2 



[(n-1) (fyWQ 2 - (n + 1) Wd*W] 



(34) 



Similarly, T can be written as 

~n— 3 



(n - l)(n - 2) (c^H/) 4 - (n + l)(n - 4)W (fyW) 2 <9jH/ 



+(n + i) 2 h/ 2 {fytrajw - (dlw) 2 } 



(35) 



Since VF^] is a polynomial function of ^ and d^W is positive, we can determine the sign 
of T' and T ' by examining the sign of polynomial in the numerators. Thus, the stability 
problem has been reduced to an algebraic one. 

Substituting the general form of W[ip) into T , we obtain 

KM 



T (r) = r 



n-2 



(36) 



where 



K\ 



m=2 
k 

m,j=2 



'2m-2 



a m (n - 1) {(3 - m)n - (m + 1)} I [ (n - p) \ ^ m_1 



p=2 



2m-2 ~\ (2j-2 

a m aj(n - 1) <j JJ (n - p) \ I JJ (n - p) 

p=2 J I p=l 



X 



j(n — 1) — (m — l)(n + 1) 



j+m-2 



(37) 



The factor other than in (1361) are manifestly positive, so the sign of K[ip] determines 
that of T . However, from fl37|) . it is clear that sign of K[ip] depends on dimensions and 
Lovelock coefficients a m . Therefore, if Lovelock black holes have the ghost instability depends 
on dimensions and Lovelock coefficients. 

However, as we will show later, T" has a negative region in even- dimensions if Lovelock 
black holes are sufficiently small. Therefore, even if T' are always positive and consequently 
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Lovelock black holes have no ghost instability, they have the dynamical instability as long 
as the ADM mass is sufficiently small. 

Substituting the explicit form of W[ip] into the formula (13^]) . we get 

T" = r^ 3 3 . (38) 

Here, the lowest and the leading term of L[ip] is 



L[ip] = ( n -l)(n-2) + --- 

4 (2k-2 ^ 

+ % | n ( n - p) 4 j ( n - ^ - - ^ - i ))^ 4fc_4 • ( 39 ) 

We note that the highest order k = [(D — 1)/2] is related to dimensions as n = 2k — 1 in odd- 
dimensions and n = 2k in even-dimensions. In odd-dimensions, the leading term disappears. 
Hence, we cannot say anything in general. Hence, we consider only even- dimensions. 

Let us examine the sign of L[ip] (tp > 0). If n = 2k, the coefficient of the lowest term 
is positive and that of the leading one of ( )39|) is negative. Therefore, L[ip] > near ip = 
and L[ip] < for large ip. This means that there exists roots of L[ip] =0 because L[i(j] is 
a continuous function. Let ipo be the lowest positive root. If ipH < ipo, then L[ip] > for 
< ip < ipHi and hence we conclude T" > for r > r#. While, Hipa > ipo, then there exists 
a region L[^] < in the range V'o < ^ < V'ff- Thus, there exists a region T" < outside 
the horizon r > r#. Therefore, black holes are stable if < and unstable if ipH > ipo- 
Since ^ becomes larger as /i becomes smaller, we conclude that there exist a critical mass 
below which black holes become unstable. 

To conclude this section, by considering tensor perturbations, we can say that 'When 
the ADM mass is sufficiently small, Lovelock black holes in even- dimensions have the ghost 
instability or the dynamical instability; that is, small Lovelock black holes are unstable in 
even-dimensions ". 

IV. STABILITY ANALYSIS FOR VECTOR PERTURBATIONS 

In this section, we consider the stability of Lovelock black holes under vector perturba- 
tions. 
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Master equation for vector perturbation is given by [281 ]; 

- d 2 ^ + V v (r)V = w 2 ^ (40) 

where 

V. (r) = -i-/*U-^) + (^ T -«)^. (41) 



Here, \& is the master variable and 7„ are eigenvalues of vector harmonics with / -/ v = t{l + 
n — 1) — 1 > 1) for k = 1 and non-negative real numbers for k = 0,-1. This equation 
is obtained provided that tensor perturbations have no ghost instability. Moreover, T{r) 
is always positive, because we assumed the conditions ©, the positivity of fi , and the 
asymptotically flat spherical symmetric branch. 

In this section, we show that Lovelock black holes are stable under vector perturbations 
as long as they have no ghost instability under tensor perturbations. In order to prove this 
statement, we again use the S-deformation approach. We define % = —d 2 /dr* 2 + V v (r) and 
the inner product (1261) . Then, as in the last section, we can find smooth function S such 
that 

(<p,H<p) = [(\D<p\ 2 + V\<p\ 2 )dr*, 
for any <p . Here, we have defined 



l 



D = ±+S, V = V v + f^-S 2 
ar ar 



Following the paper [25|, we choose S to be 



s = - f Tr^77¥ h (42) 



Then, V can be calculated as 



y = 7«-(n-l)*/r; 

n - 1 rT v ' 

In this new potential V, apparently 7^ — (n — 1)k > 0, / > and T > 0. We also assumed 

T > 0, so it is clear that V is always positive. Therefore, the inner product ((p,Htp) is 

positive for any ip. In particular, it is true for the lowest energy state, hence the lowest 

energy is positive. This implies that black holes are stable under vector perturbations. 

To summarize this section, we can say that " if there is no ghost instability in tensor 

perturbation, Lovelock black holes are stable under vector perturbations" . 
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V. STABILITY ANALYSIS FOR SCALAR PERTURBATIONS 



In this section, we examine the stability of Lovelock black holes under scalar perturba- 
tions. 



In the previous paper, we have derived the master equation for scalar perturbations 28 ]. 
Using the master variable we can write down the master equation 



8%.V + V s (r)V = tu 2 ^ . (44) 



Here, the effective potential reads 

,(rNT)' 



V s (r) = 2 ls f 



nNTr 2 

f N' 2 f t' 2 n't' 

-jjd r (fd r N) + 2/ 2 — - ^d r (fd r T) + 2f— + If — , (45) 



where we have defined 

ATI \ -Znf + 2 ls + nrf 

N(r) = = . (46) 

rvi 

For scalar perturbations, eigenvalues of scalar harmonics 7 S are given by 7 S = £(£ + n — 1) 
for k — 1 and positive real numbers for k = 0, — 1. The above master equation is obtained 
by assuming T' > 0. Hence, tensor perturbations have no ghost instability. 

Note that we will consider an asymptotically flat spherically symmetric branch with 
positive mass as in section Ull In this branch, T{r) is always positive. 

A. Criterion for Instability 

In this subsection, we show that black holes are unstable if 2T' 2 — TT" has a negative 
region outside the horizon. 

In order to prove this statement, we can use the S-deformation approach. Here, we choose 
S as 

S=fd r (\nN) + fd r (lnT) , (47) 
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then the second line of (145]) canceled and V becomes 

, (rNT)' 



27J 
2%/ 



nNTr 2 



2(7, 



UK) 



< 



nr 

Isf 

nrTT' . 



2(% - 

2T >2 



nK ) + ^npn t 2 r 



TT 



(48) 



Note that we used the assumption T > 0, T > and /i > in the last inequality. 

Now, let us prove the statement "if 2T' 2 — TT" has a negative region, black holes are 
unstable". In order to do that, we use the inequality 



2 . (<P,1l<p) 



(49) 



where Uq is the lowest eigenvalue. This inequality is true for arbitrary test function ip, so we 
choose ip as the smooth function that has compact support in the region where 2T' 2 - TT" 
is negative. Then, Uq can be bounded as 



< (ip,Hip)/(<p,<p) 
1 



< 



(<P,<P) 
1 

(<P,<P) 



dr* 



\Dip\ 2 + V\<p\ a 



dr* \D(p\' 



Is 



dr* 



f 



nrTT' 



2T 



TT 



m 



(50) 



We assume T > 0, T' > and choose as the smooth function which has compact support 
in the region 2T' 2 — TT" < 0, so the first term in (150]) must be positive and the second term 
in ( 150]) must be negative. Therefore, by taking the limit 7 S = £(£ + n — 1) — > 00, the last 
line of Eq.( l50p becomes negative, which means the lowest eigenvalue is negative. Hence, 
black holes are dynamically unstable. 

In summary, we can say that "If T has a negative region, black holes have the ghost 
instability under tensor perturbations. Even if T' is always positive, Lovelock black holes 
have the dynamical instability under scalar perturbations if 2T' 2 — TT" has a negative 
region" . Note that we can not say that black holes are stable even if 2T' 2 — TT" is always 
positive. 
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B. Instability of Small Lovelock Black Holes in Odd-dimensions 



Now let us check the sign of T' and 2T' 2 — TT" . We have already shown that Lovelock 
black holes are unstable in even-dimensions under tensor perturbations, so we concentrate 
on odd- dimensions. In order to examine the sign of these functions, it is convenient to 
express these functions by ip. The formula (|34|) reads 

T'(r)=r- 2 / [ t ] 1 , (51) 

where 

K[i/>] = {n- \){d^W) 2 - (n + V)Wd%W . (52) 
Similarly, the formula ( l35l) can be written as 

r 2n-4 

2T ' 2 - TT " = Wvr m " h (53) 

where 



M[i)] = n{n-l){d^,W) A -3n(n + l)WdlW(d^W) 



+3(n + l) 2 W 2 (dlWf - (n + l) 2 W 2 d^WdlW . (54) 



Apparently, the signs of T' and 2T' 2 - TT" are determined by K[tp] and M^]. So we have 
to check the signs of K[ip] and M[^]. 

First, we consider the sign of T' which has been already expressed by ip in fl36l) and fl37j) . 
We here consider odd- dimensions n = 2k — 1, then we have 

r , «t ^(k — 2) — 2akCXk-2(k — 1) , A , . 

Kty] = 2(k - 1) + • • • + 2 fc - 11 (fc l 1)(fc _ 2 ) V*" 4 • (55) 

where a m = a m ^Ylp=i 2 ( n ~ P)\- Furthermore, for n — 2k — 1, M[-0] can be expressed by 
ip as 

M M = 2(2* -!)(*- 1) + ... - ^"gy-V • (56) 



The most important point is that the signs of the coefficients of the leading term in ( )55l) 
and ( 1561) are opposite. 

First, we assume a^_ 1 (/c — 2) — 2akUk-2{k — 1) is negative. In this case, K[ip] is positive 
near ip = and negative for large ■0, so K[ip] has positive roots. Let the lowest root be 
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ipo- Note that ip moves in the range (0,^#]. Then, if ipn > ipo, has a negative region 
and T' does so, which means Lovelock black holes have the ghost instability. Note that tpn 
becomes larger as /i becomes smaller. Then, if a^_ 1 (/c — 2) — 2a k ak-2(k — 1) is negative, 
small Lovelock black holes have the ghost instability. 

Conversely, we assume a^_ 1 (A; — 2)— 2ak<yk-2{k — 1) is positive. In this case, the coefficient 
of the leading term of K[ip] is positive, so whether K[ip] has positive root or not is obscure. If 
this has a positive root, we can say that small Lovelock black holes have the ghost instability 
by using the same logic as the last paragraph. However, even if K[ip] has no positive root, 
M[tp] must have positive roots; this is because M[ip] becomes positive near ip = and must 
be negative for large ip. Let the lowest positive root be ipi. Using this definition, we can say 
Lovelock black holes have the dynamical instability if > ^i? namely, if Lovelock black 
holes are sufficiently small. 

In any case, assuming black holes are sufficiently small, Lovelock black hole have the 
ghost instability or the dynamical instability in odd-dimensions. 

To conclude this section, we can say that "small Lovelock black holes are unstable in 
odd-dimensions" . 



VI. EXAMPLE: EINSTEIN-GAUSS-BONNET THEORY 



In this section, to illustrate our general statement, we take Einstein- Gauss-Bonnet theory 
as an example. The Einstein-Gauss-Bonnet theory corresponds to k = 2 and generic for 
n — 3 and n — 4. 

The action for Einstein-Gauss-Bonnet theory is given by 



L = R + ^ (R 2 - + R^RxD 



(57) 



The stability of this model has been already analyzed in 23| and 25| . However, the deriva- 
tion presented here is more explicit and hence transparent. 
Substituting the concrete form of W[ip] into T', we can get 



K\ 



(n - 1) + 2a 2 (n - l)(n - 2){n - + a\{n - l) 2 (n - 2) 2 (n - 3)ip 2 . (58) 



As we can see from (1581) . T is always positive both for n = 3 and n = 4 because we are 
considering cases ip > and a 2 > 0. Similarly, substituting the explicit form of W[ift] into 
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T , we obtain 

L[ip) = (n - lf{n - 2f(n - 3)(n - 5)a 4 ^ 4 + 4(n - l) 3 (n - 2) 3 (n - 3)(n - 5)a^ 3 
+2(n - l) 2 (n - 2) 2 (5n 2 - 25n + 42)a 2 ^ 2 

+12(n - l)(n - 2)(n 2 - 3n + 4)a 2 ^ + 4(n - l)(n - 2) . (59) 

Finally, M [■?/;] which determines the sign of 2T' 2 — TT" can be calculated as 

M[<0] = -{n- l) 3 (n-2) 4 (n-3)a^ 4 + (n - l) 4 (n - 2) 3 (n - 3)a 3 ^ 3 
3 

+-(n - l) 2 (n - 2) 2 (n 2 - 5n + 2)a 2 ^ 2 

+n(ri - l)(n - 2)(n - 7)a 2 tp + n(n - 1) . (60) 

The functions L[ijj] and M[ip] give criterions for the instability of tensor and scalar pertur- 
bations, respectively. 

Let us first consider 6-dimensional Einstein- Gauss-Bonnet black holes corresponding to 
n — 4. Then, substituting n = 4 into ( I60j) . we obtain 

M[-0] = 12( — 1 + 3a 2 ^ + 9a^ 2 ) 2 > . (61) 

Therefore, black holes are stable under scalar perturbations. Next, we need to check the 
stability of Lovelock black holes under tensor perturbations. In the case of n — 4, the 
formula ( )59l) becomes 

L[ip} = -1296a 4 ^ 4 - 864a 3 ^ 3 + 1584a 2 ^ 2 + 576a 2 ^ + 24 . (62) 

Clearly, the coefficient of the leading term of L[ip] is negative. Thus, for large i/j, L[ifj] is 
negative. And, L[ip] is positive near ijj ~ because L[0] = 24 > 0. Hence, there must be a 
root somewhere between. Indeed, L[ip] becomes zero at 

^o = ^(- 1 + v / 15 + v / 10) . (63) 
6a 2 

Therefore, L[ip] is always positive in the range < ip < ipo and always negative in the range 
ip > ipo. As we explained, ip moves in the range < ip < i/jh- Therefore, if ipn > ipo, 
there exists the region T" < 0, which means Lovelock black holes are dynamically unstable. 
Furthermore, from ( Tl5|) . the inequality ipn > ipo yields 



3V6(l + v^5 + y^0) ^ 3/2 

(-i + Vi5 + VTo) 3/ 



H < — , a 2 /2 = n c . (64) 
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This proves that Lovelock black holes with the mass less than \i c are unstable in 6 dimensions. 

Now, let us consider 5-dimensional Einstein-Gauss-Bonnet black holes corresponding to 
n = 3. In the case of n = 3, the formula (15§|) reads 

L[ip] = 96a 2 ^ 2 + 96a 2 ip + 8 . (65) 

Since all the coefficient of L[ip] are positive and we are considering positive ip, L[ip] is always 
positive and hence T" is always positive. Thus, Lovelock black holes in 5-dimensions are 
stable under tensor perturbations. However, for n = 3, M[%fj\ becomes 

M[-0] = 6(1 - 4a 2 tp - Aa 2 2 i) 2 ) . (66) 

From this equation, it is easy to see that M[ifj] =0 has a positive solution; that is ipo = v ^ t ~ 1 . 
Then, M[ip] is positive in the range < ip < i/jq and negative in the range ip > i/jq. Since tp 
moves in the range < ip < i^h, 2T' 2 — TT" has a negative region if tpn > ipo- Using the 
solution for ip H (TT5]) . we can rewrite the inequality ipn > as 

(V2 + l) 2 a 2 . (67) 

Note that the lower bound came from the condition for the existence of the horizon as we 
have explained in SecJTTJ Hence, 5-dimensional Lovelock black holes with the mass in the 
above range are dynamically unstable under scalar perturbations. 



VII. CONCLUSION 



We have studied the stability of static black holes in Lovelock theory which is a natural 
higher dimensional generalization of Einstein theory. We have shown that there exists the 
instability of Lovelock black holes with small mass under tensor perturbations in even- 
dimensions and under scalar perturbations in odd-dimensions. Lovelock black holes are 
stable under vector perturbations as long as they do not have ghost instability under tensor 
perturbations. Remarkably, the instability is stronger on short distance scales, which is 
different from the usual instability for which the onset of the instability becomes a bifurcation 
point. Hence, the instability we have discussed in this paper is catastrophic in the sense that 
there is no smooth descendant. Curiously, the similar instability also appears in the Gauss- 
Bonnet cosmology 32] . In spite of this unusual nature of the instability, it is interesting 
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to investigate the fate of the catastrophic instability. This issue is very important because 
black holes lose their mass due to the Hawking radiation and eventually become unstable. 

It is worth examining the more profound meaning of this instability; that is, why Lovelock 
black holes are unstable. Especially, it is very interesting to find the reason why black 
holes have the instability under tensor perturbations in even-dimensions and under scalar 
perturbations in odd-dimensions. 

Related to the above, it is intriguing to find the thermodynamical meaning of the universal 
function T(r). As was shown in this paper, this function governs the dynamical stability 
of black holes. Therefore, if T(r) has thermodynamical meaning, the relation between 



thermodynamical 2l|, |33[ and dynamical instability might be revealed. 



It is interesting to investigate if the instability we found also exists for asymptotically 
AdS cases from the point of view of the AdS / CFT correspondence 34H36| , in particular, in 



relation to stability of holographic superconductors 
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